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Abstract 

The results of extensive Monte Carlo simulations of magnetic-field induced 

transitions in the xy model on a stacked triangular lattice with antiferromag- 

netic intraplane and ferromagnetic interplane interactions are discussed. A 

low-field transition from the paramagnetic to a 3-state (Potts) phase is found 

to be very weakly first order with behavior suggesting tricriticality at zero 

field. In addition to clarifying some long-standing ambiguity concerning the 

nature of this Potts-like transition, the present work also serves to further 

our understanding of the critical behavior at TV, about which there has been 

much controversy. 
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The possibility of novel critical behavior associated with geometrically frustrated antifer- 
romagnets has given rise to a wide variety of speculation In the cases of Heisenberg and 
xy models on the stacked triangular lattice, Kawamura has argued, by means of symmetry 
analysis and a 4 — e renormalization group expansion, in favor of a new chiral universality 
class with unusual critical exponents as determined by Monte Carlo simulations. A field 
theoretic 2 + e expansion by Azaria, Delamotte and Jolicoeur ||, however, has inspired the 
suggestion that such systems exhibit nonuniversality where first-order, mean-field tricritical 
or 0(4) criticality can occur depending on unspecified details of the model (also see Ref. f§|). 
Two very recent Monte Carlo simulations on the Heisenberg model have yielded different 
conclusions dependent upon details of the analysis. Bhattacharya et al || maintain that this 
system exhibits 0(4) universality whereas our own work supports the idea of a very weak 
first-order transition with possibly a pseudo-critical region |7[ giving 0(4) critical exponents. 
Two similar numerical studies of the Ising model, although yielding substantially the same 
results, have led to the speculation of yet another new universality class by one group 0, in 
contrast with our intrepretation |5| that the previous suggestion of standard xy universality 
]T0| is confirmed. It is also of interest to note that the Monte Carlo results which led to the 
proposal of a new universality class associated with the frustrated pyrochlore antiferromag- 
net Jill have been re- interpreted in support of a first-order transition ||12j| . It is becoming 



increasingly clear that the results of direct numerical simulations of frustrated spin systems 
can be difficult to interpret. The present work represents an attempt to reveal the critical 
behavior of the xy model on a stacked triangular lattice by application of a magnetic field H. 
The analysis of these Monte-Carlo simulation results is guided by expected behavior based 
on symmetry arguments of a less controversial nature and lend support to the proposal of 
tricriticality (or an extremely weak first-order transition) for this system. 

This work was inspired by the extensive study of Lee et al |13[ who examined the xy 
antiferromagnet on a triangular lattice (un-stacked) in an applied magnetic field. At H=0, 
the transition exhibits Kosterlitz-Thouless behavior but the field breaks rotational symmetry 
and transitions involving true long-range spin order occur. For H > 0, but not too large, a 



colinear phase is stabilized with the symmetry of the 3-state Potts model. The real 3-state 
Potts model exhibits a continuous transition in two dimensions (2D) and, after numerous 
numerical simulations over the past twenty years, it appears to be generally accepted that the 
transition is weakly first order for the 3D lattice |14|,|15j. For the so-called continuous 3-state 
Potts model, an effective Landau-Ginsburg- Wilson (LGW) Hamiltonian is constructed which 
contains a term cubic in the order parameter [Jig]. Within mean-field theory, such models 



yield a first-order transition, independent of space dimensionality. The transition in 2D is 
thus driven to be continuous by critical fluctuations, with known critical exponents verified 
by the work of Lee et al. Some analyses of renormalization-group and series expansions for 
the 3D case indicate a transition which may be continuous but most studies favor the first- 



order scenario [14,17]. The conclusion of an earlier Landau-type analysis, which partially 



included effects of fluctuations, is that models of this type may exhibit either a continuous 



or first-order transition depending on relative parameter values |Tj|. Although contrary 
results would have been surprising, it was not a priori certain that the transition to the 
3-state ordered phase considered here would be first order. We find convincing evidence 
that the transition is indeed weakly first order. This study represents the first Monte Carlo 
simulation of a model equivalent to the continuous 3-state Potts Hamiltonian in 3D. 
We study the Hamiltonian 

n = h £ Si • Si (i) 

<ij> 

where the spins lie in the basal plane, Jy is the interplane interaction, Jj_ > indicates 
the antiferromagnetic coupling which is frustrated for the triangular geometry, < i,j > and 
< k,l > represent near-neighbor sums along the hexagonal c axis and in the basal plane, 
respectively, and the field is applied in the basal plane direction x. The magnetic order 
realized by this model can be conveniently described in terms of a spin density expressed as 
a low-order Fourier expansion M 



sir 



m + Se iQr + S*e- iQ - r (2) 
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where m is the uniform component induced by the magnetic field, Q is the wave vector, and 
the complex polarization vector can be written in terms of real vectors, S = S a + zS&. The 
120° spin structure known to occur at zero field below the Neel temperature Tjv is described 
by a period-3 basal-plane modulation, along with a helical polarization for S. At H=0, the 
critical behavior is independent of the sign of the interaction Jii. For Jn < there is no 
interplane modulation but for Jm > a period-2 structure is stabilized. This difference gives 
rise to a term cubic in S in a Landau-type free energy, or LGW Hamiltonian, only for the 
case J\\ < 0: F 3 ~ (m • S)S • S) + c.c. For a linearly polarized spin density, one can write 
S = S r e 4<?i , where S r ||m is real, to get a contribution F 3 ~ mSfcos(3(p). The three (Potts) 
states arise from the three inequivalent choices of the phase angle = nir/3, where n is an 
integer. Higher-order terms in the free energy can stabilize <fi = (2n + 1)tc/6 depending on 



field and temperature values ||19|| . In this case the cubic term is zero. 

In order to determine the magnetic-field temperature phase diagram, standard Monte 
Carlo simulations were performed on the Hamiltonian (1) with Jii = — 1 and J± = 1 for 
lattices Lx L x L with L=12-24. Runs of 1 - 2 x 10 4 Monte Carlo steps (MCS) per spin 
were made, with the initial 4—8 x 10 3 MCS discarded for thermalization. Boundary lines were 
estimated as in our previous work [[2tJ . Not surprisingly, the result shown in Fig. 1 is similar 
to the 2D case studied by Lee et al. In particular, in addition to the paramagnetic phase 
1, there are two ordered phases, 6 and 9, with a linear polarization of the spin density and 
an elliptical phase 7 (phases are numbered following a previous convention ^lf). At H=0, 
the expected 120° spin configuration was observed, with the Neel temperature T/v — 1.45 
in agreement with Kawamura @ and Ref. [pO[j . These phases have the same symmetry as 



determined in the 2D case. Phase 6 with S||H is the 3-state Potts phase discussed by Lee 
et al. 

A molecular field treatment of this model P^j , which yields results independent of di- 
mension, gives a phase diagram where the linear state 6 is absent and two critical lines 
which merge at T/y. This is somewhat puzzling since the cubic term responsible for the 
stability of phase 6 occurs in the expansion of the free energy, in powers of s(r), which can 



be derived within this mean-field theory. Such an expansion gives an expression for F(m,S) 
which is identical in structure to what is obtained for a phenomenological Landau-type free 
energy (also mean-field) |Jl9"| , p0f . The only difference is that in the latter case each of the 
(six) fourth-order terms have an independent coefficient Bi (since each term is an indepen- 
dent invariant), whereas in the molecular-field treatment all these coefficients are equal. We 
find that by making just one of the Bi different from the others, a phase diagram with the 
correct structure (Fig. 1) is found, where the 1-6 boundary is first order as expected fl23 



Molecular- field theory appears to somewhat accidentaly exclude phase 6. This model also 
yields the result that the 1-9 phase boundary represents a line of continuous transitions since 
the phase angle <fi approaches the value ir/6 at this boundary line so the cubic term is not 
relevant. 

The criticality of the 1-6 transition boundary was studied at two points using the 
Ferrenberg-Swendsen histogram method of analysing Monte Carlo data. This technique 
is well-suited for the study of transitions which may be very weakly first order, particularly 
when used to determine the internal-energy cumulant [ |TT| , p5[| 

U(T) = 1 - - < E A > I < E 2 > 2 . (3) 
3 

This quantity exhibits a minima near Tjy, which achieves the value U* = | in the limit 
L — > oo for continuous phase transitions. In the case of a first-order transition, U* < | 
is expected. The histogram method may also be used to determine precisely the location 
of extrema near Tjy which occur in other thermodynamic functions. These are expected to 
demonstrate simple asymptotic volume dependence in the case of first-order transitions or 
have L-dependence governed by critical-exponent ratios in the case of continuous transitions. 
In addition to U, results are given here for the specific heat (C), staggered susceptibility (%) 
as well as the logarithmic derivative of the order parameter (V) which is equivalent to pB] 

V(T) =<ME>/<M>-<E>, (4) 

where the relevant order parameter M is defined as in Ref. [|20| . Simulations were performed 
at H=0.7 and H=1.5 on lattices with L=12-33. Thermodynamic averages were made using 



1 x 10 6 MCS for the smaller lattices and up to 2.6 x 10 6 MCS for the larger lattices, after 
discarding the initial 2 — 3 x 10 5 MCS for thermalization. If necessary, several runs at 
different T were made to ensure that the extrema of the desired function occured close to 
at least one simulation temperature. The estimated critical temperatures are 1.488(2) and 
1.522(2) for H=0.7 and 1.5, respectively. 

Asymptotic scaling of the extrema with volume for the case H=0.7 demonstrated in Figs. 
2-4 is consistent with a first order transition 0. An indication that it is only weakly first 
order is revealed by noting that the estimate U* = 0.66660(3) from Fig. 2 is more than 
an order of magnitude closer to | than the value determined for the 5-state Potts model 



in 2D pv |, considered to be one of the weakest first-order transitions known. It is also of 



interest to note that the value U* = 0.6460(2) was determined by Fukugita et al [|15] for the 



discrete 3-state Potts model on a cubic lattice, also known to be only weakly first order p8 
Similar first-order bevavior was found in the data at H=1.5. Evidence that the 1-6 transition 
becomes more strongly first order as H increases is given by the estimate U* = 0.66643(3) at 
this higher field value. Histogram data were also taken at H=0.7 for the 6-7 transition, which 
has a critical temperature estimated to be 1.420(2). In a manner similar to our analysis of a 
continuous transition on the stacked triangular lattice ||, scaling consistent with the Ising 
universality class is evident. This is the expected result as this transition involves only the 
order parameter S&y. Details of these results, along with histogram analyses of the other two 
transition lines, will be presented elsewhere ||23|| . Preliminary results for the 1-9 transition 
indicate this transition is continuous, as it is within the phenomenological Landau-type 
model discussed above. 

The corresponding phase diagram for the case of antiferromagnetic interplane coupling 
has only one linear state with S_LH and two critical lines emanating from T/v, which are 
transitions of xy (S\) and Ising (Z 2 ) universality p0]| . This structure nicely reveals the 
Z2 x S\ symmetry of the order parameter for the transition at T/v 0. The results of the 
present study suggest that this picture does not occur in the case of J\\ < (or for Jy > 
with an applied field staggered along the c axis). 
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The Monte-Carlo histogram simulations of this work give clear indication that the present 
version of the continuous 3-state Potts model exhibits a very weak first-order transition in 
3D. In addition, our results suggest that this transition becomes more weakly first order as 
the field is lowered. T/v thus appears to have characteristics in common with a tricritical 
point, a possibility suggested by Azaria et al ]3[]. This scenario is made somewhat ambiguous 
by the fact that T/v is also a multicritical point where more than one phase meet. The con- 
ventional Monte Carlo simulations of Kawamura were used to estimate critical exponents 
associated with his proposed chiral universality class which are not very different from those 
expected for mean-field tricriticality (a possibility that was considered in Ref. @). In view 
of these results and the recent ambiguity in interpreting Monte Carlo data for the frustrated 
Heisenberg model |5|||], as well as the discussion by Peczak and Landau of pseudo-critical 
behavior associated with weakly first-order transitions, our results are consistent with the 
transition at TV for the xy model being tricritical or very weakly first order. It appears that 
only histogram Monte-Carlo simulations at T/v which are very extensive (long runs on large 
lattices) have the possibility to add new information on this problem. 
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FIGURES 

FIG. 1. Phase diagram determined by standard Monte Carlo simulations (points with error 
bars). Indicated are the paramagnetic phase 1, phases 6 and 9 having colinear order and phase 
7 with an elliptical (chiral) spin structure. Squares at H=0.7 and 1.5 indicate boundary points 
determined by highly accurate histogram analyses. Solid and broken lines are guides to the eye 
and indicate first and second-order transitions, respectively. 

FIG. 2. Scaling of the energy-cumulant minima with volume for the 1-6 transition at H=0.7 
The straight line represents a fit to the four largest lattice sizes. 

FIG. 3. Scaling behavior of the specific heat maxima with volume as in Fig. 2. 

FIG. 4. Scaling behavior of the maxima of the susceptibility \ an d logarithimic derivative of 
the order parameter V as in Fig. 2. 
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